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ABSTRACT 

We characterize spinning dust emission in the warm ionized medium by comparing templates of 
Galactic dust and Ha with the 5-year maps from the Wilkinson Microwave Anisotropy Probe. The 
Ha-correlated microwave emission deviates from the thermal bremsstrahlung (free-free) spectrum 
expected for ionized gas, exhibiting an additional broad bump peaked at ^ 40 GHz which provides 
^ 20% of the peak intensity. We confirm that the bump is consistent with a modified Draine & 
Lazarian (1998) spinning dust model, though the peak frequency of the emission is somewhat lower 
than the 50 GHz previously claimed. This frequency shift results from systematic errors in the large- 
scale modes of the 3-year WMAP data which have been corrected in the 5-year data release. We show 
that the bump is not the result of errors in the Ha template by analyzing regions of high free-free 
intensity, where the WMAP K-band map may be used as the free-free template. We rule out a pure 
free-free spectrum for the Ha-correlated emission at high confidence: ~ 27(7 for the nearly full-sky fit, 
even after marginalizing over the CMB cross-correlation bias. We also extend the previous analysis 
by searching the parameter space of the Draine & Lazarian model but letting the amplitude float. 
The best fit for reasonable values of the characteristic electric dipole moment and density requires an 
amplitude factor of ~ 0.3. This suggests that small PAHs in the warm ionized medium are depleted 
by a factor of ~ 3. 

Subject headings: diffuse radiation — dust, extinction — ISM: clouds — radiation mechanisms: non- 
thermal — radio continuum: ISM 



1. INTRODUCTION 

Following formu l ation of a physical model 
(jPraine fc LazarianI Il998bl . hereafter DL98), spin- 
ning dust has been invoked to explain the dust- 
correlated microwave emission seen by the Wilkinson 

Micro wave Anisotrop y Probe (WMAP) from 23-41 

GHz (de Oliveira-Costa et all 120041: iFinkbeined l2004 
Finkbeiner, Langston. & Minter [20041 : iBoughn fc Poben 
2007.: ,Dobler fc Finkbeiner .2008bl ). Given the WMAP 
data at 94 GHz, the dust-correlated emission at 23-41 
GHz was anomalously high to be explained by a purely 
t hermal emission mech anism. 

iBennett et al.l (|2003D attempted to obviate the incon- 
sistency by suggesting that the first year data from 
WMAP did not require a spinning dust component to 
construct an internally consistent foregrounds model, so 
long as the dust-correlated emission was actually dust- 
correlated synchrotron at low frequencies. However, 
when the 3-year WMAP data (23-94 GHz) were com- 
bined with lower frequency data from Green Bank (8,14 
GHz; iFinkbein er. Langston. fc Minter '2004), Tenerife 
(10,15 GHz; Be Olivcira-Costa ct al. 2004), and other 
telescopes (19 GHz: IBoughn fc Pobeil 120071) . it became 
clear that the dust-correlated emission at ~ 20 GHz did 
not behave spectrally like synchrotron emission, and was 
instead consistent with a spinning dust spectrum. 

The principal source of ambiguity in interpreting the 
low frequency WMAP data as evidence for spinning dust 
was that the spectrum of the dust-correlated emission 
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falls from 23-41 GHz. The implication is that either 
the peak in the spectrum is near or below the lowest 
WMAP band or that the emission is synchrotron, and it 
is only with data at lower frequency that the spinning 
dust spectrum is recovered. In the DL98 spinning dust 
models however, the peak frequency of the spinning dust 
emission varies depending on model parameters (e.g., the 
grain size, geometry, ambient density, etc.), and at higher 
temperatures with smaller grain size is expected to be 
^ 40 GHz. This regime is termed the "warm ionized 
medium" (WIM, T - 8000 K) in DL98 compared to 
the "cold neutral medium" (CNM, T 100 K) emission 
that peaks at ~ 20-30 GHz and is thought to be the 
source of the dust-correlated emission at low frequencies 
in WMAP. 

Dobler fc Finkbeiner (2008b, hereafter DF08b) showed 
that the ambiguity is eliminated and that the turnover 
in the spinning dust spectrum is recoverable within the 
WMAP frequency range using Ha-correlated emission as 
a tracer of WIM spinning dust emission. Typically, Ha 
is used as a tracer of free-free (thermal bremsstrahlung) 
emission from ~ lO^'K gas in regions where extinction 
from dust is minimal. However, because a map of Ha rep- 
resents an emission measure (density squared integrated 
along the line of sight) and because spinning dust emis- 
sion is generated via dust grains that are spun up by 
collisions with ions (which should also be roughly pro- 
portional to density squared integrated along the line of 
sight), the Ha emission should also spatially trace WIM 
spinning dust emission. In the analysis of DF08b, the 
peak frequency of the emission was found to be ~ 50 
GHz, which was slightly higher than the ~ 40 GHz pre- 
dicted by DL98. 

In this paper, we show that the evidence for spinning 
dust in the Ha-correlated emission has strengthened with 
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the release of the 5-year data, due to both decreased noise 
as well as corrections made to large-scale fluctuations 
that were present in the 3-year data. 

2. TEMPLATE FITS 

Microwave emission mechanisms in the WMAP fre- 
quency range can be separated into four broad categories: 
free-free, thermal dust, spinning dust, and synchrotron 
emission.* Free-free emission is generated by electron-ion 
collisions which produce thermal bremsstrahlung in ion- 
ized gas; thermal dust refers to electric dipole emission 
from thermal fluctuations (e.g., lattice vibrations) in the 
electric charge distribution in the grain; and spinning 
dust is electric dipole radiation from the smallest dust 
grains which are excited into rotational modes through a 
variety of coUisional mechanisms. Synchrotron emission 
consists of a soft component (T oc v^^^ with (3$ ~ 3) 
originating from supernova shock accelerated electrons 
which spiral around the Galactic magnetic field and a 
hard component (/? « 2.5) centered on the Galactic cen- 
ter and extending roughly 20 degrees. This hard com- 
ponent has been termed the "ha z e" and its origin re- 
mains uncertain (jFinkbeineil 12004 iDobler fc Finkbeh^ 
HoOSa). 

Each of these emission mechanisms is approximately 
traced by maps of the sky at other frequencies, described 
below. This external information about the spatial struc- 
ture of the foregrounds makes possible a multi-linear fit 
of the spectrum of each foreground. This multi-linear 
regression may be performed over the whole sky, or in 
selected regions to study the spectral variation of each 
component. We now briefly describe the templates used , 
and the fitting method (see lDobler &: Finkbeineil [200881 
for more details). 

2.1. The templates 

For fitting purposes, we distinguish between the 
soft synchrotron component originating from supernova- 
shock accelerated electrons, and a harder synchrotron 
component in the inner Galaxy, possibly with a different 
physica l origin. Soft s y nchro tron emission is well traced 
by the iHaslam et all (|1982| ) 408 MHz map. Though 
the spectral index of this emission is e xpected to vary 
slightl y from place to place in the sky, iLa Porta et al.l 
(|2008f) showed that from 408 to 1420 MHz, T cx v'f^^ 

with f3s ~ 3; A similar spectral index was derived by 

iPage et al.l (|2007D using polarization in the WMAP data. 

The WMAP data, even in the first year, contained 
evidence of a harder synchrotron component in the in- 
ner Galaxy, which has become known as the "haze" be- 
cause of its relatively featureless morphology (jFinkbeineil 
120041) . Dobler & Finkbeiner (2008a, hereafter DF08a) 
interpret the haze as a separate physical component in 



* IDraine & : LazariajJ l|1999l) point out that if interstellar dust 
contains an appreciable fraction of ferrimagnetic or ferromagnetic 
materials, thermal fluctuations in the magnetization would gener- 
ate appreciable amounts of magnetic dipole emission at frequen- 
cies V < lOOGHz, which would be a fifth type of CMB foreground. 
However ICassasus et al.l pOOSl ) find that this type of emission is 
inconsistent with the morphology of cm emission from the p Oph 
cloud, implying that magnetic materials are not abundant in inter- 
stellar dust. Thus, we do not consider it further in this paper; we 
will use the term "thermal" emission to refer specifically to emis- 
sion from thermal fluctuations in the electric dipole moments of 
grains. 



the template fit described in DF08a and below, and the 
spectrum of this emission was found to be T cx v~^^'^ 
with Ph ~ 2.5. The template we use for this compo- 
nent is T cx l/rgai where rgai is the distance from the 
galactic center in degrees. Though systematics related 
to contamination of any CMB template by residual fore- 
grounds leads to uncertainty in the derived soft and hard 
synchrotron spectra, we note that neither the presence 
nor the spectrum of the haze has any significant effect 
on the conclusions in this paper. For a detailed study 
of microwave synchrotron in the Galaxy see Dobler & 
Finkbeiner (2008c, in preparation). 

Our thermal dust template is the Finkbeiner, Davis, 
& Schlegel (1999, hereafter FDS99) two-component 
model of thermal dust, evaluated at 94 GHz. The 
FDS99 map uses column densi ties an d tempera- 
tures from IRA S /ISS A (jW heelock eTan [r994l and 
COBE/DIRBE (|DIRBE Exp. SuppTfT99g) given by 
Schlegel, Finkbeiner, & Davis (1998; SFD), constrained 
to fit COBE/FI RAS data in the microwave and sub-mm 
(jFIRAS Exp. Su pp. 1997). The preferred FDS99 model 
(model 8) gives 1^ cx i^^ '^ B^{Tdnst)- Converting this to 
antenna temperature, T cx with (3d k, 1.55 in the 
WMAP bands for typical Tdust values. We use the same 
template for the diffuse, cold neutral medium (CNM) 
spinning dust which has a peak frequency ~ 20 GHz. 

For a free-free and WIM spin ning dust t e mplat e we 
use the Ha map assembled by iFinkbein e? (I200 3D us- 
ing d ata from thr ee surveys: VTSS (jPcnnison et al.l 
'1998), SHASSA (Gau stad et aD l200ll l. and WHAM 
yiaffncr et al. 2003). The template is corrected for dust 
extinction assuming the dust and ionized gas are uni- 
formly mixed along the line of sight. This approxima- 
tion fails in regions of very high dust column density, so 
we mask out regions in our fit where the SFD extinction 
at Ha is ^(Ha) = 2.QhE{B - V) > 1 mag. Addition- 
ally, we mask out all point sources listed in the WMAP5 
point source list as well as the LMC, SMC, M31, Orion- 
Barnard's Loop, NGC 5090, and the HII region around 
C Oph. This mask covers 22.2% of the sky and is shown 
in Figure [2] below. Roughly 1/4 of the masked pixels are 
point sources. 

2.2. CMB Estimators 

Finally, we also need an estimator for the CMB since it 
contributes a large variance at WMAP frequencies (and 
in fact constitutes the biggest source of noise when de- 
riving foreground properties). We use the six CMB esti- 
mators described in DF08a and summariz e their features 
in Ta ble [H As noted in DF08a (see also iHinshaw et ahl 
l200l . every CMB estimator is contaminated to vary- 
ing degrees by foregrounds which leads to a systematic 
bias in the inferred foreground spectra. In the case of 
"internal linear combination" (ILC) type estimators, in 
which the WMAP maps in each band h are weighted and 
summed with weight coefficients Cfc that are chosen to 
approximately cancel the foregrounds while preserving 
unity response to the CMB, the bias is proportional to 
the chance spatial correlation of the true foregrounds (as 
opposed to the templates) with the CMB (see DF08a). 

The spectral cross correlation coefficients given in e.g. 
Fig. m depend on the CMB estimator used, but because 
we fit out a CMB spectrum in the interpretation of the 
correlation spectra (see f|5|), the final results for the Ha- 
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TABLE 1 



CMB type 


description 




ILC coefficients 








K 


Ka 


Q 


V 


W 


1 


published WMAP Syr ILC 






N/A 






2 


ILC using WMAP Syr Kp2 coefficients 


0.134 


-0.646 


-0.377 


2.294 


-0.40S 


3 


minimum variance ILC over unmasked, fit pixels 


0.048 


-0.603 


0.342 


0.S30 


0.683 


4 


WMAP W band minus (thermal dust + free-free model) 






N/A 






5 


minimum variance ILC with thermal dust model presubtracted 


0.182 


-0.638 


0.044 


0.32S 


1.086 


6 


TOH Syr map* 






N/A 







Note. — The different types of CMB estimators used in the fits. Because of contamination of these estimators by fore- 
grounds, the inferred foreground spectra vary depending on the estimator used. Our preferred estimator is CMBS, which 
consists of an ILC with a minimum variance (over unmasked pixels) weighting of the thermal dust model presubtracted 
WMAP data. One of the attractive features of this map is that it has well understood noise properties. 

"available at http://lambda.gsfc.nasa.gov/ 
''Tegmark et al. priv communication 



correlated spinning dust do not depend on our choice of 
CMB estimator. 

2.3. Template fit procedure 

To infer the spectra of the individual foreground com- 
ponents, we use the multi-hnear regression template fit 
described in DF08a. We construct the template ma- 
trix P, whose columns consist of the templates described 
above, and derive the spectrum of the template corre- 
lated emission, a^, where a* denotes the spectrum of fore- 
ground i and a.b is the vector of correlation coefficients 
for band b. For each band, we solve the matrix equation, 
Pat = Wf,, (1) 
for ab, where Wf, is the WMAP map for band b. To deter- 
mine the best fit eh, in Equation [1] we evaluate the "pseu- 
doinverse" P+ (see DF08a) and note that af, — P+Wf, 

minimizes the quantity = \ Psib — Wf,|^ where the 
sum is over unmasked pixels. Divid ing P and Wh by a, 
(the mean noise in e ach band^ , as in 'Bennett et"al]|200 
iHinshaw et al.|[2007[ ). this solution for a^ minimizes, 

2 



\Pab - Wf,| 



xl 



(2) 



The total — J2b xl for the fit is given by the sum over 
the 5 WMAP bands. We also define the residual maps 
rb = W6 - Pat. 

With our template fits, the spectrum of each fore- 
ground emission component a* is left completely uncon- 
strained, though it is assumed to be constant across the 
sky. To appropriately account for the CMB, we also force 
g^CMB = units of thermodynamic temperature. 

3. FOREGROUND CROSS-CORRELATION SPECTRA 

Our foreground spectrum fits are characterized by the 
cross-correlation coefficient a as well as the formal er- 
ror in the fit in each band. These errors take into ac- 
count both the noise in each WMAP band and morpho- 
logical correlations between the different templates. By 
our definition of a above, the derived (cross-correlation) 
spectra are in I^, units of kJy/sr per template unit — 
Rayleighs for Ha-correlated emission and mK for FDS99- 
correlated emission. Both the data and templates are 
mean-subtracted over unmasked pixels to remove sensi- 
tivity to zero-point levels in the maps. 

^ The WMAP scan strategy visits some pixels, e.g. the ecliptic 
poles, far more than average. Because we are in the systematics- 
dominated limit, we use the average noise in each band to avoid 
over-weighting certain parts of the sky. 



3.1. Ha and Dust Correlated Emission 

Figure [1] shows our derived cross-correlation spectra 
for the Ha- and FDS99-correlated microwave emission 
from 23 to 94 GHz for our six CMB estimators. The 
dust-correlated spectrum exhibits the familiar thermal 
tail from 94 to 61 GHz and then the rise from ^^50 to 
23 GHz from anomalous emission. While the statistical 
error bars on the spectra are very small due to both the 
high sensitivity of WMAP as well as the large sky cov- 
erage (- 150,000 pixels for HEALPix N^idc = 128), the 
systematic effects from the contamination of the CMB 
estimators by foregrounds are significant. 

As identified in DF08a and DFOSb, the spectrum 
of Ha-correlated emission does not follow the free-free 
power law as expected. Instead, there is a bump in the 
spectrum with a peak frequency i/p « 40 GHz. This 
bump is present for all CMB estimators. We note that 
it cannot be generated by a contamination of the CMB 
estimator by foregrounds (since this bias has the spec- 
trum of the CMB) and DFOSb argued that this bump is 
most easily explained by a WIM spinning dust compo- 
nent that is traced by the Ha map. We address both of 
these considerations in detail below. 

3.2. Changes from 3-year to 5-year data 

There are two notable differences when comparing our 
results using the 3-year versus 5-year WMAP data. First, 
the inferred peak frequency Vp w 40 GHz is significantly 
lower than the 50 GHz seen in the 3-year WMAP data 
(see DF08b). This shift in peak frequency is primarily 
due to spurious large-scale power that was present in the 
3-year data but has been corrected in the 5-year data re- 
lease. Figured] shows that the amplitude of these fluctua- 
tions can be quite large. In particular, the V band dipole 
is clearly apparent and has peak values from ~ — 7/iK to 
- 10/iK. 

As described above, we have mean subtracted the maps 
over unmasked pixels to perform the template fits, and 
so our results are not sensitive to zero point offsets. How- 
ever if there is spurious low £ powe r, as in the 3-year data , 
our fit results will be affected. IHinshaw et al.l ()2008f ) 
point out that these large-scale fluctuations in the 3-year 
data were due to imperfect characterization of the instru- 
ment gain. While the amplitudes of the fluctuations are 
not large enough to have a significant impact on CMB 
analyses, Ha-correlated emission mechanisms (free-free 
and WIM spinning dust) are sufficiently subdominant at 
WMAP frequencies that ^ 10/iK offsets become impor- 
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Fig . 1.- — Cross-correlation spectra for the Ho {left), and FDS99 dust (right) templates using various estimators of the CMB anisotropy 
(see i|2.2l l. Although the formal error in the fit is small, in both cases the unknown degree of CMB contamination introduces a highly 
covariant uncertainty. In each case, the spectra are the same up to addition of a CMB spectrum of unknown amplitude. 



tant. Because the 61 GHz band (V band) was the most 
affected, and because the WIM spinning dust peak; fre- 
quency is in the range 30-50 GHz, it is not surprising that 
our best fit peak; frequency has shifted. Had the spurious 
low-^ power been morphologically similar in all bands, 
the peak frequency would have been similar between the 
3- and 5-year data. The problem was exacerbated by 
the GMB estimators, which necessarily contained some 
complicated combination of these fluctuations. 

The second difference between the 3- and 5-year anal- 
yses is that, for a given GMB estimator, the ILC coeffi- 
cients (defined in Equation |4] below) have changed. Al- 
though this does not affect our interpretation of the fit 
results (see Q , the spectrum of Ha-correlated emission 
has also changed by more than the formal fit uncertain- 
ties. Since this is due primarily to changes in the co- 
efficients used when constructing a GMB estimator, the 
extent of the contamination of these estimators by fore- 
grounds has changed from year three to year five. For 
example, in the case of a simple ILG which minimizes 
the variance over unmasked pixels (GMB3) with respect 
to the ILG coefficients Cb, 



0, 



where 



TlLC — 



Cb^Vb 



(3) 



(4) 



is the ILG map in thermodynamic mK and the ILG co- 
efficients are constrained to sum to unity to preserve 
response to the GMB. These coefficients are sensitive not 
only to the large-scale fluctuations mentioned above, but 
also to the noise in the measurements themselves since 
Wb — c -I- ff, -I- n;,, where c, ff,, and n;, are the GMB, 
total foregrounds in each band, and noise in each band 
respectively. Thus the signal variance in the noise also 
contributes to the determination of Ct,. This is illustrated 
by the limit of large measurement noise, in which case 



the minimum variance linear combination would have 

3.3. Ha as a tracer of spinning dust 

For the 3-year WMAP data, DF08b presented an ar- 
gument for why the Ha map should trace spinning dust 
emission at WMAP frequencies. Briefly, the Ha intensity 
is proportional to the emission measure (EM), J n^di = 
EM, where is the electron density and the integral is 
along the line of sight. We will show below that eval- 
uating the DL98 spectrum with parameters appropriate 
for the WIM, the total emission is indeed roughly pro- 
portional to rig. That is, including other (de-)excitation 
mechanisms such as plasma drag, far-infrared emission 
from the grains, etc., the behavior of the total emission 
with density implies that the Ha EM map should trace a 
component of the spinning dust emission in the Galaxy. 

4. SPINNING DUST MODEL 

In this section, we summarize the spinning dust model 
used in our analysis of the 5-year Ha- and dust-correlated 
spectra below. The model for excitation mechanisms 
and dipole moment distributions is described in detail 
in DL98. We find that the most relevant parameters are 
the density, nn and the typical electric dipole moment of 
1 nm grains, /ip. 

We use th e upda ted grain size distribution given in 
iDraine fc Lil (|200l . Our grain size distribution is a 
double lognormal distribution characterized by a mean 
and width for the small and large grain populations — 
g) = 4 and 20 A and aa = 0.4 and 0.55 respectively 
Draine fc Lil[200l . This size distribution has a lower 
cutoff of ainin = 3.55A which corresponds to iVatoms = 20 
Garbon atoms. The dipole moment distribution for the 
grains consists of three delta functions so that 50% of 
the grains have /i — /io(a/l nm)'^/^, 25% have /i half as 
large, and 25% have /i twice as large. 

The spinning dust emission in the DL98 model comes 
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Fig. 2. — The difference between the WMAP 5-year and 3-year data in the 5 bands. Left: WMAP5— WMAP3. Updated beams and gains 
lead to differences in the plane, however large scale gradients (and in some cases dipoles) are seen in each band. Center: WMAP5— WMAP3 
masked with the pixels that are used in the multi-linear regression fit (see J2]l. Right: WMAP5— WMAP3 smoothed to 30° FWHM and 
masked. Here the large-scale fiuctuations are clearly evident. In particular, the peaks range from ^ —8 /.tK to ~ 15 fiK. All maps were 
mean subtracted over unmasked pixels prior to subtraction and smoothing. 



TABLE 2 



Environment 


nH(cm ^) 


fio (debye) 


Tgas (K) 


Tdust (K) 




xm 


X 


y 


WIM 


0.01-0.6 


0.5-12 


3000 


20 


0.99 


0.001 


1.0 


0.0 


CNM 


5.0-40.0 


0.5-12 


100 


20 


0.0012 


0.0003 


1.0 


0.0 



Note. — Spinning dust parameters for the "warm ionized medium" (WIM) and "cold 
neutral medium" (CNM) environments. For a detailed definition of the parameters shown here 
see DL98. 



from polycyclic aromatic hydrocarbons (PAHs) of very 
small size, A^atom ~ 100 or less. Small asymmetries in 
the grain geometry result in non-zero dipole moments 
that are taken to scale with the number of atoms as 
M oc V A'atom • As these grains are spun up by various ex- 
citations processes (see DL98), they emit electric dipole 
radiation with total power P (x uj'^, where uj/2tt is the 



angular rotation frequency and is tens of GHz.^ 

The parameters which characterize the environment 
and grain properties, with values appropriate for the 
WIM and CNM (cold neutral medium), are summarized 

^ An improved treatment of the spinning dust problem by Ali- 
Hamoud & Hirata has made a few refinements of the DL98 model 
but those modifications do not significantly affect the results of this 
paper (Ali-Hamoud & Hirata priv. comm.) 
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Fig. 3. — Upper left: Spinning dust emission per H« intensity for 
WIM parameters, keeping fia constant and varying njj- Increasing 
rill increases the peak frequency but changes the total power per 
Ha little. Upper right: Same, but for constant njj and varying fio- 
Decreasing /iq also increases the peak frequency but decreases the 
total power per Ha. Note that in this part of (nH,A*o) parameter 
space, the shape of the spectrum is nearly unchanged, but shifts 
in frequency and amplitude. Middle panels: Contours of the peak 
frequency (left) and peak amplitude (right) in the (riHiMo) plane. 
Lower panels: Same, but for CNM parameters. See the discussion 
in [|4]for details. 



in Table [21 The top two panels of Figure [3] shows ex- 
ample spinning dust spectra for WIM parameters while 
holding the density nu fixed while varying the charac- 
teristic dipole moment /xq and vice-versa. These spectra 
are in units of Iu/Ihoi, or total intensity per intensity 
at Ha with Tgas = 3000 K (our justification for using 
this gas temperature, which is lower than the commonly 
used 8000 K, is given in ijS.ip . For fixed as nn is 
increased, the spectrum shifts up in frequency with very 
little change in overall power per Ha. Since I Ha oc 
in fully ionized environments, the total spinning dust in- 
tensity goes roughly like . This point will be explored 
in more detail shortly, but it illustrates that an EM map 
like the Ha map should trace spinning dust emission. 

For fixed tt-h, as /ip is increased Figure [3] shows that the 
spectrum increases in amplitude and decreases in peak 
frequency. This can be understood in terms of the radi- 
ated power per grain which goes as /i^ : grains with larger 
dipole moments radiate away power faster and are thus 
harder to spin up leading to a lower peak frequency, but 



they also emit more total power leading to an increased 
amplitude. 

The lower four panels of Figure [3] are contour plots 
showing the effects of varying ripj and /.tg on the am- 
plitude and peak frequency of the spinning dust spec- 
trum. Contours are shown for both the WIM and CNM 
conditions. The contours of peak frequency show that, 
if the amount of dust per H is assumed unknown (i.e., 
the amplitude is allowed to vary), then there is a strong 
degeneracy in the peak frequency between njj and /iq- 
The implication of this strong degeneracy is that a mea- 
surement of the spinning dust spectrum with sparse fre- 
quency sampling that only gives information about the 
peak frequency (as is the case with the WMAP data be- 
low, see Figure [7]) cannot uniquely constrain these two 
parameters. On the other hand, if the grain abundance 
is assumed known, than the degeneracy is broken and 
the parameters are more tightly constrained. 

Figure [4] shows the behavior with nu of the emissivity 
per EM for both the total emission, J judv, as well as vji, 
evaluated at the five WMAP frequencies. The parame- 
ters of the WIM and CNM environments are also given 
in Table [2] while those for the reflection nebula (RN) 
and photo-dissociation region (PDR) are shown on the 
plot. In the limit of low rtn, all four environments exhibit 
roughly ji, cx njj behavior. In this limit, the dominant 
excitation and de-excitation mechanism are absorption 
and emission of infrared photons and thus the total emis- 
sion scales linearly with density. This is in contrast to 
the argument presented in DF08b who argued that in the 
lowest density regions, the grains are in the episodic limit 
and so the total emission should scale as n^. That argu- 
ment neglected the importance of photon interactions. 

While the linear nn scaling persists through all densi- 
ties shown here for the RN and PDR regions, the WIM 
and CNM both flatten to ji, oc behavior at nn ~ 0.01 
and 5 cm~'^ respectively. In this regime, the dominant 
spin up mechanism is collisions with ions and the 
scaling persists for four orders of magnitude in njj for 
the WIM (~ 10~^ — 10^ cm~^). As the density increases 
further in the CNM, the spectrum turns over and again 
goes roughly like j^, cx nn- At these high densities, both 
spin up and damping are dominated by thermal processes 
(collisions with ions and atoms). However, it is interest- 
ing to note that over the density ranges shown in Table[2l 
both the WIM and CNM models scale as n\ indicating 
that a map that scales with EM like the Ha map should 
trace WIM and CNM spinning dust emission. 

5. THREE COMPONENT FOREGROUND SPECTRA 

Given that the Ha map should trace both the free- 
free and a spinning dust emission component, and given 
the fact that our foreground spectra are necessarily con- 
taminated by a component with the spectrum of the 
CMB (flat in thermodynamic AT) from the CMB cross- 
correlation bias described in DFOSa and S[2l we interpret 
the Ha-correlated emission in Figure [T] as a three com- 
ponent spectrum. As in DFOSb, we fit a free-free plus 
CMB plus DL98 WIM spinning dust model. 
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Fig. 4. — Spinning dust emission per for the WIM (upper left) CNM (upper right), RN (lower left), and PDR (lower right). Eacli 
panel contains vjv/ri^ for the 5 WMAP bands: 23 GHz (solid), 33 GHz (dashed), 41 GHz (dash dot), 61 GHz (dash triple-dot), and 94 GHz 
(long dash). The upper soUd Une is the integrated emissivity, J judu/rv^. Some panels include a line slope of -1 (dotted), corresponding 
to jv/riYi constant with njj- This limit is generally approached when photons dominate, either because of low density or intense radiation 
field. The solid vertical line marks the value chosen by DL98 for each environment. 



to the measured cross-correlation spectrum I^. Here, 
pic is the "planck correction" factor which converts ther- 
modynamic AT to antenn a temperature at frequency v 
(|FIRAS Exp. Supp.1 [19971 ). We choose the parameters 
njj and /xg, generate a spinning dust spectrum, and min- 
imize = '^iilv — , where Ui are the errors 
in the fit in each band i, over the parameters Fq, Cq, 
and Dq. We allow Dq to float in the fit because of uncer- 
tainty about the PAH size distribution and abundance 
in the WIM. We concentrate on CMB5, but our results 
are minimally changed for the other CMB estimators, it 
is only the value of Cq which is significantly affected. 

5.1. Ha- Correlated Emission 



Figure [5] shows the results of our three component fit 
for the Ha-correlated emission. From the figure, it is 
clear that a spinning dust spectrum with nn and /ig set to 
0.15 cm~^ and 3.5 debye respectively (and with the WIM 
parameters in Table ^ fits the data remarkably well. 
This value for /io is comparable to the earlier estimates in 
DL98. This spectrum has a peak frequency of Vp = 37.2 
GHz and Z?o fit coefficient of 0.34. It is interesting to note 
that the value of Dq is less than unity. This may indicate 
either that PAHs are depleted in the diffuse WIM by a 
factor of ~ 3 or that the grain size distribution differs 
from that in the CNM from which the values of the log- 
normal mean and width of our grain size distribution 
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Fig. 5. — Upper: The Ha-correlated spectrum from 23 to 94 
GHz (points) using CMB5 and the three-component fit to the data 
(dashed line). The three components, free-free, spinning dust, and 
CMB, are shown separately. The amplitude of the free- free compo- 
nent yields Tgas ~ 3000 K, while the amplitude of the spinning dust 
component is less than unity indicating that either PAHs are de- 
pleted in the WIM or the grain size distribution is altered. Lower: 
A zoom in of the spectrum showing that, although the bump is 
subtle, the statistical significance is very high. 

were derived (Drai ne fc Lill2007t ). 

From the amplitude of the Fq fit coefficient, we can 
derive a gas temperature in the WIM which we find to 
be Tgas ~ 3000 K. This is the origin of the choice of Tgas 
in the WIM spinning dust model. This value for Tgas is 
much lower than the temperature inferred from emission 
line ratios (e.g., see Madsen et al. 2006). We stress that 
this low gas temperature is not the result of assigning 
some of the intensity to a spinning dust component. Ig- 
noring the bump and merely fitting the data with a free- 
free only spectrum, still gives a gas temperatur e below 
5000K . Using a technique similar to ours, Davie s~et al.l 
(|2006f ) also report anomalously low gas temperatures 
based on the free-free to Ha ratio in the WMAP 1-year 
data and DF08b found similar temperatures in the 3-year 
data. The origin of this discrepancy is unclear, though 
it may result from temperature variations i n the WIM 
along the line of sight as discussed in lHeilesI (|2001h .^ We 
will explore implications of a two-component tempera- 
ture model in future work. Lastly, we find that the to- 
tal intensity in Ha-correlated spinning dust emission is 
roughly 25% of that in free-free emission at 41 GHz. 

The upper panel in Figure [5] shows that the bump in 
the spectrum is indeed subtle, but the zoom in beneath 
reveals that it is highly statistically significant. Since 
we minimize over three amplitudes {Fq, Dq, and Co), 
set the values of nn and /xq by hand, and fit to five data 
points, the number of degrees of freedom is formally zero. 
However, Figures [6] and [7] show that the rises quickly 
as we deviate from these parameters. In the left panel 

^ Another possibility is that our template fit of the spectrum 
is dominated by the strongest Ha features which naturally have a 
lower temperature. 



of Figure [6l Ax^ contours are shown in the (Do,Co) 
plane and it is clear that the null hypothesis that the 
Ha-correlated emission is simply a linear combination of 
free-free plus a CMB bias (i.e., Dq = 0), is ruled out at 
very high confidence. In fact we find that the Dq ~ 
case is ruled out at ^ 27a. The right panel shows 
contours in the (£'o,^'p) plane as Vp is varied by hand. 
This parameter is less well constrained but is broadly 
consistent with i/p 35-40 GHz. 

Figure [7] shows x^ and Dq contours in the (nH,/io) 
plane for the full-sky fit. The contours show that 
there is a long narrow valley in the surface reflecting 
the degeneracy in peak frequency shifts resulting from 
varying nn and /zq (see Figure Since varying /iq af- 
fects the overall amplitude of the spinning dust compo- 
nent in units of Iv/ iHa while nu does not (because of the 
density squared behavior of the spinning dust emission 
at these rin), the value of Dq changes as we move along 
the narrow valley of the x^ surface. It is important to 
note that small deviations from these parameters quickly 
lead to large x^ values. Further, given this set of environ- 
mental conditions and grain sizes, Dq = 1 is ruled out; 
again, implying either PAH depletion or modified grain 
size distributions in the WIM. 

5.2. The Gum Nebula 

In addition to our nearly full sky template fit, we have 
also carried out smaller regional fits. Of particular in- 
terest are the four regions of the Gum Nebula shown in 
Figure 4 of DF08b. This is a region that is very bright 
in Ha but has little in the way of thermal dust emission 
or synchrotron. The null hypothesis is that the spec- 
trum of this region should correspond almost entirely to 
a free-free spectrum. 

Figure [8] shows the three component model fit as well 
as the contour plots for the four regions of the Gum Neb- 
ula. The bump in the Ha-correlated spectrum persists 
for all four regions, though the spinning dust amplitudes 
Dq and CMB correlation bias amplitudes Co vary slightly 
from region to region. There is some evidence that the 
peak frequency Vp also varies slightly from region to re- 
gion, but it is the least well constrained parameter and 
all four regions indicate roughly Vp ~ 35-45 GHz. From 
the Ax^ contours in the (Dq, Cq) plane, we find that the 
null hypothesis of Dq = is ruled out at high confidence 
in all regions. 

Due to the larger error bars in the template fit (the 
smaller regions include fewer pixels), the x^ valley in the 
('T'HiMo) plane is somewhat broader than in the full sky 
fit. Nevertheless, all of the spectra are consistent with 
a WIM spinning dust spectrum with rtn ~ 0.1 cm^"^, 
fiQ ^ 2 debye, and Dq ~ 0.4. 

5.3. Thermal Dust-Correlated Emission 

As shown in Figure [1] dust-correlated emission falls 
from 23 GHz to 60 GHz and it was this departure from a 
thermal spectrum which initially lead to speculation that 
this "anomalous" dust component originated from spin- 
ning grains. In this section, we fit our three component 
model to this data but replace the WIM spinning dust 
model in Equation [5] with a CNM spinning dust model 
and the free-free component with a thermal dust spec- 
trum T = To X (i//94 GHz)'3o. We choose Pd = 1-7 but 
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Fig. 6. — Left: Ax^ contours in the {Do,Co) plane. The null hypothesis that the Ha-correlated emission follows a free-free power law 
with a contamination by the CMB spectrum is ruled out to very high significance (x^ ~ 740 for 3 degrees of freedom at Dq = 0). Right: 
contours in the (foi'^pcak) plane. Though the peak frequency is not constrained to high accuracy, i^peak ~ 37 GHz is roughly the best 

good fit to the data (x^ = 50.5), especially considering 
that there are few free parameters. However, the right 
hand panel of Figure [9] indicates that, although the am- 
plitude is approximately correct (To w 1.0), for the CNM 
parameters there is no point in parameter space which 
yields a significantly better fit. 

The poor fit is due to the fact that the 23, 
33, and 41 GHz data points exhibit roughly pow- 




erlaw behavior (T 



-2.83 



) SO that the spinning 



Fig. 7. — The (dashed) and Dq (solid) contours in the 
{nil, fio) plane. There is a deep narrow valley in the x^ surface 
representing the near degeneracy in effects on the spinning dust 
peak frequency model with increasing njj and decreasing fj,o- The 
values nil = 0.15cm~^ and = 3.5 debye produce an acceptable 
fit with Do ~ 0.3. Although Dq does vary along the narrow valley 
in x^, Do = 1 does not yield an acceptable fit, indicating that either 
the PAHs are depleted in the WIM or the grain size distribution is 
altered. 

note that our results are not significantly changed for the 
range 1.5 < /3_d < 2. 

The left hand panel in Figure [5] shows the results 
of fitting a CNM spinning model to the WMAP dust- 
correlated spectrum. The spinning dust model has njj = 
11 cm^^ and /iq — 5.5 debye. From the figure, it is clear 
that this particular spinning dust spectrum is not a very 



dust spectrum has too much curvature to fit the 
data. This powerlaw behavior has lead several au- 
thors to misidentify this emission as "dust-correlated 
synchrotron" , howe ver both WMAP polarization data 
(|Kogut et al.l I2007D and data at lower frequencies 
over large areas jFinkbeiner. Langston. fc Minted 120041 : 
Ide Oliveira-Costa et al.ll2004f) argue stronglv against the 
synchrotron hypothesis. Furthermore, we note that, 
along any line of sight, there may very well be regions 
with different environment and grain properties leading 
to a superposition of spinning dust spectra. This degen- 
eracy between synchrotron and a superposition of spin- 
ning dust spectra can only be broken with sufhcient fre- 
quency coverage. In particular full-sky, high-resolution 
maps between 5 and 15 GHz could eliminate the ambi- 
guity 

Lastly, we point out that the anomalous dust- 
correlated emission in the CNM could originate from 
an emission mechanism that is neither synchrotron nor 
spinning dust. For example, iDraine fc LazarianI ()1999f ) 
suggest that emission from magnetized dust grain could 
contribute to the total emission in the lower frequency 
WMAP bands. However, in the case of the WIM, the 
fact that the Ha map and the DL98 spinning dust mod- 
els both scale roughly as density squared suggests that 
the bump in the Ha-correlated emission originates from 
spinning dust and not D L99 magnetic dust. Further- 
more, the observations of iCassasus et al.l ()2008[) appear 
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Fi g. 8. — The same as Figures [5] [B] and [7l except for the four regions of the Gum Nebula defined and illustrated in lDobler fc Finkbeined 
(|2008b') . The null hypothesis that the Ha-correlated emission is well represented by a free- free spectrum with a CMB contamination is 
ruled out to high significance in all regions (x^ > 14 for 3 degrees of freedom at Do = 0). 



to strongly limit the contribution of magnetic grain ma- 
terials to the microwave emission from dust. 



6. 



ALTERNATIVE EXPLANATIONS FOR THE 
Ha-CORRELATED BUMP 



In this section we will address alternative explanations 
for the origin of the bump in the Ha-correlated emis- 
sion. The alternatives listed here are not meant to be 
an exhaustive list of possibilities but rather constitute 
potentially important features of the maps, our fitting 
procedure, etc., and thus require careful attention about 
whether they can artificially generate a bump in our spec- 
trum. For each of these, we argue why the bump is more 
likely to be explained by spinning dust. 

6.1. Bandpass Effects 



Although the spectra presented in previous sections 
consisted of five data points, the individual channels of 
WMAP have broad {Ai'/i/ « 0.2) bandpasses that must 
be taken into account. For a given emission mechanism, 
the total intensity in a WMAP channel is the spectrum 
of that emission integrated over the bandpass of that 
channel. 



TtOt 



J ijj{v)dv 



where is the total intensity in band h and io{v) is the 
bandpass for band h normalized so that J uj{i')dv = 1. 



Thus, for a free- free spectrum (/^ 



the ob- 



served spectrum will be slightly modified from a straight 
powerlaw when plotted against the CMB weighted band 
centers (see below). We find that the deviations are at 
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Fig. 9. — Left: A three component model fit to the FDS99-correlated emission spectrum from 23 to 94 GHz. The fit includes a thermal 
dust component (T oc i'^ '^), a CNM spinning dust component (rifj = llcm~^ and fig = 5.5 debye; see Table [2J, and a CMB component. 
Close inspection reveals that the curvature in this single component spinning dust spectrum is incompatible with the data and over-predicts 
the amplitude at 61 GHz. Right: The same as Figure [7] except for CNM parameters and fitting to FDS99-correlated emission. Although 



the amplitude is roughly Dq = 1, the values are unacceptably high. In fact, for the model in the left hand panel, 
reasonable CNM parameters, a single component DL98 spinning dust model is ruled out at high confidence. 



: 50.5. Thus, for 



the level of ~ 1%. Nevertheless, this total intensity, eval- 
uated at each band, is what must be used to fit the de- 
rived cross-correlation spectrum. 

Additionally, the WMAP data themselves must be con- 
verted from thermodynamic AT to Ii, in each band as- 
suming a CMB weighted band center,* 



(7) 



where the spectrum of the CMB /^^^ oc v'^plc{v) (i.e., 
flat in thermodynamic AT) . All of the fits presented in 
previous sections take the bandpass weighting into ac- 
count, and we find that the changes to the results are 
also at the ^ 1% level. 

6.2. Dust Extinction and Masking Effects 

It has long been known that, although a map of Ha 
is an excellent tracer of free-free emission from ionized 
gas, it is not perfect. In particular, extinction of Ha by 
dust makes the map a poor tracer of free- free i n regions of 
high d ust column density. To account for this lFinkbeineil 
()2003l ) (whose map we have used in this paper) apply a 
dust correction w hich assumes un iform mixing between 
the dust and gas. [Dickinson et al.l {2003} take a slightly 
different approach by matching slices through the galac- 
tic plane to an extinction model plus dust column density 
model and attempting to directly measure the extinction. 

The effects of dust on the Ha map could manifest 
themselves in two ways. First, recall that our mask is 
generated by avoiding regions where the extinction due 
to dust at Ha is E{B — V) > 1 mag. It is difficult to 

* By design the WMAP receivers measure energy and so the 
natural data unit is f I^du, however since they are then calibrated 
off of the CMB dipole, the data are presented in thermodynamic 
AT. 
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Fig. 10. — Foreground spectra derived via our multi-linear 
regression template fit as a function of the extinction cut used 
for the mask. The bump in the Ha-correlated emission persists 
for all mask definitions indicating that it is not due to improperly 
corrected extinction in the Ha map. 



envision that our mask could have a sufficiently large 
impact on our derived spectrum to produce the bump, 
but it is instructive to consider the possibility. In Fig- 
ure [TU] we show the derived Ha-correlated spectrum for 
a wide range of dust cuts (i.e., for numerous masks). As 
can be seen in the figure, the bump feature is robust to 
variations in the mask. 
The second potential issue is that perhaps the Ha map 
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has been over corrected for dust so that it is actuaUy a 
rather poor tracer of free-free, potentially leading to the 
bump as a spurious artifact. To test this hypothesis, we 
performed our template fit with an Ha map that was not 
corrected for dust extinction at all. Even in this extreme 
case of wnrfer-correction, we find that the bump persists. 
We conclude that the bump does not result from errors 
in the extinction correction of the Ha map. 

6.3. CMB Cross-Correlation Bias 

The CMB cross-correlation bias described above (and 
more extensively in DF08a) is both large and ubiquitous 
in CMB foreground analyses. Because the effect on the 
derived spectrum of each component can be substantial, 
it is important to rule out the possibility that this bias 
is producing the bump. This possibility represents the 
"null hypothesis" discussed in H5.ll The basic question 
is, can the derived cross-correlation spectrum be fit by 
a simple linear combination of free-free and CMB spec- 
tra? As we showed in that section, this Dq = case is 
ruled out at very high confidence with a P^r degree 
of freedom of 245 (for 3 degrees of freedom at Dq = 0) . 

6.4. Cross- Correlation Between the Templates 

If there is any chance spatial cross-correlation between 
the templates, than this will impact the shape of the de- 
rived spectra. Furthermore, if one of the templates is 
contaminated by Ha morphology or vice-versa, than the 
spectrum of each foreground will be "contaminated" by 
the others. Again, a clear example of this would be poor 
dust correction of the Ha map adding dust-like morphol- 
ogy to that map. Another example would be reflection 
of Ha photons created in the plane off of dust at high 
latitudes, which would also imprint a dust-like morphol- 
ogy on the Ha map. There is also the fact that the WIM 
emits thermal radiation as well and so there is some Ha 
morphology present in the FDS99 map. 

DFOSb showed that the chance spatial cross- 
correlations between the templates used in the analysis 
is sufficiently small that there should be relatively little 
contamination of one foreground spectrum onto another. 
We explore this point in greater detail in Appendix |^ 
where we specifically address the fact that WIM features 
appear at a low level in the dust map. 

6.5. Variations of Tgas with Position 

Gas temperatures in the interstellar medium are known 
to vary with position. For example, there is evidence that 
there is a decrease i n Tgas towards the Galactic center 
(|Quireza et al.l[2006h . However, even though relatively 
small variations in Tgas lead to large changes in the Ha 
to free-free ratfo, th e effect is simply an amplitude shift 
(lValls-Gabaudl[T99l . Thus, whe n we fit the nearly full 
sky, we are averaging over many regions with slightly 
different amplitudes, but all with spectra cx i/^'^-^^. Thus, 
our resultant spectrum should be cx indicating 
that a superposition of many regions with varying Tgas 
would not produce a bump in the average spectrum. 

6.6. Model Independent Identification of the Anomalous 

Emission 

Rather than attempting to address every possible con- 
tamination mechanism that could result in a bump in 



the Ha-correlated spectrum, it is possible to identify the 
bump without the use of any template fitting and only 
using the WMAP data. The top two panels of Figure 
[11] show the WMAP data at 23 GHz (with CMB5 sub- 
tracted) as well as the Ha map. We identify two regions 
that are bright in Ha and also bright at K band (23 
GHz): C, Oph near (£,6) = (6,23) degrees and the Gum 
Nebula near £ = 260 and -22 < 6 < 17 degrees . If the 
Ha map traces only free-free emission, then subtracting 
K band scaled by a free-free spectrum from Ka, Q, V, 
and W band should completely eliminate these features. 

The bottom grid in Figure fTTl shows K band in the left 
most column and then the result of subtracting K band 
scaled by a free-free spectrum from Ka, Q, V, and W 
bands (columns 2, 3, 4, and 5 respectively) for both the 
Gum Nebula region (row 1) and the C Oph region (row 
2). Since the spectrum of free-free is harder than syn- 
chrotron and softer than thermal dust, synchrotron emis- 
sion is highly over-subtracted and thermal dust emission 
is under-subtracted (note the prominence of thermal dust 
emission in column 5). What we find is that with a free- 
free spectrum, there is a clear excess in the southern 
Gum Nebula region and particularly in C, Oph. The im- 
plication is that, C Oph and the Gum Nebula contain 
an emission component which is harder than free-free 
from 23 GHz to 41 GHz, but softer than free-free from 
41 GHz to 94 GHz. Furthermore, this spectrum is not 
soft enough at 23-61 GHz to be synchrotron, yet is vastly 
brighter than the expected thermal dust emission, so it 
cannot be any of the standard foreground components. 
Rather, it has a spinning dust type spectrum with peak 
frequency near 40 GHz. 

7. DISCUSSION 

We have repeated the multi-linear regress ion template 
fit outlined in Dob ler fc Finkbeinej ()2008al DF08a) on 
the 5-year WMAP data and found that the anomalous 
"bump" at ^ 40 GHz in the Ha-correlated emission spec- 
trum persists and has in fact become more pronounced. 
The most significant changes in the spectrum from 3-year 
to 5-year data come from large scale modes (particularly 
a large dipole at 61 GHz V band) of amplitude ^ 10 /iK 
that were corrected in the 5-year data. 

We interpret this bump as a "warm ionized medium" 
(WIM) spinning dust component that is traced by the 
Ha map. We find that for WIM parameters, the to- 
tal spinning dust emission per grain scales roughly as 
the ambient ionized gas density squared (see Figure H]). 
Therefore, a map of emission measure (EM = J nidi) 
like the Ha map would be expected to trace this compo- 
nent. 

Any derived foreground spectrum is subject to the 
CMB correlation bias (see DFOSa) from chance spatial 
correlations of the foreground emission with the CMB. 
Thus, we fit a three component model consisting of free- 
free, WIM spinning dust, and CMB spectra to the de- 
rived Ha-correlated emission in which the amplitude of 
each component is allowed to float. The best flt model 
gives a value for the ion density of the ambient medium 
of riH ~ 0.15 cm~^ and a characteristic dipole moment 
for the grains of ^ 3.5 debye referenced to 1 nm grains. 
However, there is a strong degeneracy between these two 
parameters as they have similar effects on the peak fre- 
quency. 
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Fig. 11. — Top left: WMAP K band minus CMB5. Top right: The Ha map. First column: A cut out of K band around the Gum 
Nebula (row 1) and Oph (row 2). Second through fifth column: K band subtracted from Ka, Q, V, and W bands using a free-free u^^-^" 
spectrum. Clearly there is residual emission in both the southern Gum Nebula and the Oph cloud using a free-free spectrum indicating 
the need for a separate component. This component is harder than free-free from 23 to 41 GHz and softer than free-free from 41 to 94 
GHz, which is indicative of a spinning dust spectrum. Although there are significant regions of over-subtraction (due to the substantial 
synchrotron component at low frequencies), the Ha-correlated emission is clearly under-subtracted by a free- free spectrum. 



The amplitude of the spinning dust spectrum must be 
reduced by a factor ^ 0.3 compared to the model to fit 
the data, possibly indicating that PAHs are depleted in 
the WIM. From the amplitude coefficient of the free- free 
component, we use the Ho; to free- free ratio to infer a gas 
temperature of ^ 3000 K. Near the peak frequency at 41 
GHz, the total spinning dust emission is roughly 20% of 
the free-free emission. We also find that the bump per- 
sists in smaller regional fits of the Gum Nebula as well 
with roughly the same amplitude. There is some evi- 
dence that the peak frequency varies with region, but this 
parameter is the most poorly constrained in our analysis. 

As in numerous other studies, we also find the now fa- 
miliar rise in the thermal dusi-correlated spectrum from 
41 to 23 GHz from anomalous dust-correlated emission. 
An attempt to fit a three component model of thermal 
dust, "cold neutral medium" (CNM) spinning dust, and 
CMB spectra to the dust-correlated emission yields a 
poor fit. However, we expect that, along any line of sight, 
there will be multiple regions with diflferent environmen- 
tal properties so that the low frequency dust-correlated 
emission may represent a superposition of spinning dust 
spectra. Although other possibilities exist for this emis- 
sion (magnetic dust and "dust-correlated synchrotron" 
being two), we point out that since the Ha map is a den- 
sity squared map, and since WIM spinning dust emis- 
sion goes roughly as density squared, we believe it is 
very likely that the bump in the Ha-correlated emission 
represents a spinning dust component. Magnetic dust 



emission would not correlate with Ha in this way. 

We have shown that the bump cannot be explained 
by numerous possible systematic uncertainties in our fit. 
For example, the bump cannot be the result of the CMB 
correlation bias because a fit to the spectrum with a two 
component free-free plus CMB model yields a very bad 
= 740 with three degrees of freedom. It also cannot 
be due to poor correction for dust extinction of the Ha 
map since the bump persists even if the template fit is 
performed with a completely uncorrected map, nor is it 
sensitive to our choice of mask for which we cut on dust 
extinction of Ha. Variations in gas temperature with po- 
sition and cross-correlation between the templates used 
in the fit are also incapable of producing the bump. 

Perhaps the most striking visualization that the bump 
is not an artifact of the fitting procedure is by subtracting 
the K band WMAP data scaled with a free-free spectrum 
from Ka, Q, V, and W bands. This reveals that in re- 
gions that arc very bright in K band and Ha which were 
thought to be completely free-free dominated at these 
frequencies (e.g., the Gum Nebula and ( Oph), there is 
a positive residual at 33, 41, and 61 GHz. Scaling by a 
free-free plus bmnp spectrum instead removes the resid- 
ual. The implication is that there is emission coming 
from these regions that has a spectrum that is harder 
than free-free from 23-41 GHz and then softer than free- 
free from 41-94 GHz, which resembles a spinning dust 
spectrum. 

The interpretation of the dust-correlated anomalous 
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emission as spinning dust (or magnetic dust) has been 
plagued by the fact that ancillary data sets at lower fre- 
quencies are needed to distinguish that emission from 
synchrotron. Here we have shown that a turnover in the 
spinning dust spectrum is completely recoverable within 
the WMAP frequency range. The only surprise is that 
it correlates with the Ha map. In hindsight, this is not 
very surprising since the Ha map traces density squared 
emission, though it does suggest that there is no "true" 
template for the spinning dust. Rather, it is potentially 
ubiquitous and probably will not correlate precisely with 
any map of the sky. 
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APPENDIX 

A. IMPERFECT TEMPLATE MORPHOLOGY 

Here we consider the effects of "contaminating" one of the templates with another, and determine the effects on 
the resultant spectra. We phrase this in terms of a template (denoted with primes) and a "true" map that perfectly 
traces the emission for the component we're solving for. We consider two cases. First, we consider that the Ha map is 
not a perfect tracer of free-free emission and show that this does not mix the spectrum of other foregrounds into the 
Ha-correlated spectrum. That is, the bump in the spectrum is not due to contamination effects from other foreground 
spectra. Second, we consider the case in which the dust template, FDS99, is tracing both a CNM and WIM spinning 
dust component at low frequencies. This case does have an impact on the spectrum, but we show that it cannot 
produce the observed bump and that it does not significantly impact our main conclusions. 
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Fig. 12. — An illustration of the contamination of FDS99 by Ha morphology. The left panel shows the FDS99 dust map centered on 
the Gum Nebula {{£, b} = {260, 0}, a strong Ha feature). The top middle panel shows the best fit FDS99 = egL + SqH where L is a map 
of H I emission and H is the Ha map. The top right panel shows that the difference between the model and the FDS99 map contains 
very little Ha morphology (lower right panel). Nevertheless, the effect is subtle as shown in the residual plus SqH map in the lower middle 
panel. Our fit of the level of contamination of the FDS99 map by Ha morphology is Sq ~ 4.5 X 10~^ mK/R. 
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Fig. 13. — Left: the spectrum of Ha-correlated emission using FDS99 — 5 X Ha as a CNM template for various val ues of 5. Note that 
5 Ax 10-^ mK/R is roughly required to appropriately remove Ho morphology from FDS99 as shown in Figure [121 As expe cted, as S 
is increased, the spectrum from Figure [5] is contaminated with FDS99-correlated spectrum in Figure |9] (see text, Equation IA8I I. Middle: 
the spectrum of Ho-correlated emission assuming that FDS99 exactly traces thermal dust emission with T(^^J5t = FDS99 X {u/94 GHz)'^-^. 
For 5 = mK/R, the bump in the Ho-correlated spectrum represents the difference of WIM and CNM spectra, for 5 = 4 X 10~^ mK/R, 
the spectrum is almost entirely WIM, and as S is increased further, the bump represents an addition of WIM plus increasingly more CNM 
spectra. Indeed, the peak frequency moves to the lower CNM value as 6 gets large. Right: The (dashed) and Do contours in the (njj, ^J.o) 
plane for the spectrum shown in the middle panel with (5 = 4 X 10~^ mK/R. With this small level of contamination, our results for the 
best fit riH, AtOi a^nd Dq values are almost completely unchanged (cf. , Figure [7j. 



A.l. Ha as an imperfect tracer of free- free 

For the template fits the equation we are trying to solve, neglecting the CMB and the associated cross correlation 
bias and considering only two foregrounds, is 

aiH' + a2D ^tiH + t2D, (Al) 

for ai and a2, where H is the true Ha map which we assume exactly traces free-free emission, D is the dust map, ti and 
t2 are the true amplitudes (at an arbitrary frequency), and H' is the observed Ha map (which may be contaminated 
in some way). 

By multiplying Equation lAll by H' and solving for oi, we obtain 
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and likewise by multiplying by D and solving for 02 we find 

{DH) {DH') 

a2^t2+hj^-ai^^. (A3) 

To show that the Ha-correlated emission spectrum is not affected by the spectrum of the dust-correlated emission, we 
need to show that ai 7^ 01(02, t2)- Plugging Equation IA3I into Equation I A2I we get, 

{ HH') _ {DH){DH') {DH'){DH ') 

{iP^ ^ {D^){H'^-) ^ {D^){H'- 
or 

/ {DH'){DH') \_ ({HH') {DH){DH') \ 
{D^){H'^) )^^'[{H'^) {D^){H'^) )■ 

Equation IA5I shows two things. First, oi does not depend on either 02 or t2, implying that the spectrum of Ha 
correlated emission does not depend on either the spectrum or the inferred spectrum of dust correlated emission. But, 
there is an overall multiplicative normalization, that depends on cross correlations of if, H' , and D. Note that as 
H —t H' , ai ^ ti. Equation IA5I is the general case; but now suppose we assume a specific form for if', namely we 
mix in a bit of the dust map so that H' = H + eD. This can originate from imperfect dust correction of the Ha map 
or scattering of Ha photons off of dust for example. Then the right hand side of Equation I A5I reads . 

^ {HH') {DH){DH') \^ f{{H'~eD)H') {D{H' - eD)){DH') 



(ii'2> (D2)(i/'2) J {H'^) (L>2)(iJ'2) 

{DH') {DH'){DH') {D'^){DH') 



-U 1 



' (iJ'2) {D^){H'^) " (D2)(iJ/2) 



Comparing the term in parentheses with the left hand side of Equation IA5I shows, that, in the case where H' is a 
linear superposition of H and D, ai = ti. That is, we recover the true spectrum. The bottom line is that, if the Ha 
map does not exactly trace the free-free morphology, the spectrum will be affected, but only at most by a constant 
multiplicative factor (i.e., it will not produce a bump in the spectrum). 

A. 2. FDS99 as a tracer of both WIM and CNM spinning dust 

So far we have assumed that FDS99 is a tracer of exclusively the CNM dust, as if WIM dust only appeared as 
Ha-correlated emission. While the WIM dust indeed has a hotter color temperature (e.g. in the 60/100 micron ratio) 
than CNM dust, it certainly also appears in the FDS99 map. Therefore, in our regression analysis, fitting the CNM 
component with the FDS99 template will inadvertently absorb some of the WIM emission and affect the inferred 
spectrum of the WIM. The sense of the effect is to bias the Ha-correlated WIM emission with some amount of the 
CNM spectrum, pushing the WIM spectrum to higher peak frequency and lower amplitude. In this section we estimate 
the extent of this bias. 

In order to calculate an estimate of the contamination of FDS99 by Ha morphology, we model the FDS99 map as a 
linear combination of maps of H I emission and Ha emission, 

D^eoL + SoH, (A7) 

where again, H is the Ha map, and L is the iKalberla et aTl (|2005f ) H I map. With this notation, Sq has the units of 
mK/R. We solve the above equation over unmasked pixels (where the mask consists of the point source mask described 
in §2.11 plus a masking of all pixels for which the dust extinction at Ha is 2.6 5E(B — ^) > 0-5 mag) and find that 
^0 ~ 4.5 X 10^^ mK/R. Although this value varies slightly from region to region (|Heiles et al.lll999l ). the contamination 
is very small and extremely subtle (see Figure fT2|). Nevertheless, based on our value of <5o. Ha morphology does indeed 
"contaminate" the FDS99 map since the large grains in the WIM emit thermal ra diatio n at 94 GHz as well. 

To calculate the affect on the inferred spectrum, we can use the same procedure as ^A.ll but now taking D' D+SqH. 
It is straightforward to show that the recovered Ha-correlated spectrum becomes, 

ai^ti-Sot2. (A8) 

That is, the inferred iia-correlated spectrum is contaminated by the FDS99-correlated spectrum at a level that is 
given by So- Assuming that H exactly traces both the free-free emission and the WIM spinning dust emission and 
that D exactly traces CNM emission (thermal and spinning dust), the implication is that the bump in the inferred 
spectrum of the Ha-correlated emission represents the WIM spinning dust spectrum minus some amount of the CNM 
spinning dust spectrum, determined by how much the Ha morphology has leaked into the FDS99 map via large grain 
thermal emission in the WIM. 

The left panel of Figure [13] shows the level of contamin ation of the Ha-correlated spectra by taking D' — FDS99 — Sx 
Ha as a tracer of CNM dust. As expected from Equation lASl for large values of S, the Ha map is being over-subtracted 
from FDS99 and the resultant spectrum is highly contaminated by a CNM thermal plus spinning dust spectrum (cf.. 
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Figures [5] and [9|) . For S — niK/R on the other hand, some of the WIM spinning dust bump is being absorbed in the 
/^'-correlated spectrum and so the bump in the Ha-correlate d sp ectrum represents a WIM spinning dust minus some 
CNM spinning dust spectrum. Based on the fit in Equation IA7[ the S = Sq ^ A x 10^^ mK/R case is a nearly pure 
mixture of free-free plus CMB plus WIM spinning dust spectrum. 

However, since the WIM does also emit thermal dust radiation, if we assume that FDS99 is a perfect representation 
of the thermal emission at 94 GHz, then removing the WIM (Ha) morphology from FDS99 implies that the Ha- 
correlated spectrum also contains that thermal emission. We can estimate the amplitude of this effect by assuming 
that the WIM thermal emission has a dependence on frequency with /3 = 3.7 (in intensity units, /3 = 1.7 in antenna 
temperature) and adding S x H x {v/di GHz)'^ '' back into the Ha-correlated spectrum. This is equivalent to setting. 



The resultant spectra are shown in the middle panel of Figure [T51 In this case, the bump is ideally due to only spinning 
dust emission from the WIM for S = A x 10^^ mK/R. For larger values of 6 the CNM spinning dust spectrum begins 
to leak into the fit so that by (5 = 10 x 10~^ mK/R, the spinning dust bump has moved towards the lower CNM value, 
and the amplitude has increased since the WIM and CNM spectra are now both present in the data. 

Lastly, the right panel of Figure [T51 shows and Dq contours in the (nHiMo) plane for a free- free plus CMB plus 
WIM spinning dust fit to the S — 4 x 10~^ mK/R spectra in the middle panel of Figure [Ol These contours indicate 
that the best fit values of ?ih and /io are not significantly affected by the presence of WIM morphology in the FDS99 
map for our estimate oi 5 — Ax 10~^ mK/R. Nevertheless, the exact level of PAH depletion measured by our technique 
(Do) will depend on precisely how one estimates the WIM leakage in FDS99. We defer a detailed analysis of this to 
future work. 
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